Doping driven Small-to-Large Fermi surface transition and d-wave superconductivity 

in a two-dimenional Kondo lattice 
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We study the two-dimensional Kondo lattice model with an additional Heisenberg exchange be- 
tween localized spins. In a first step we use mean-field theory with two order parameters. The first 
order parameter is a complex pairing amplitude between conduction electrons and localized spins 
which describes condensation of Kondo (or Zhang-Rice) singlets. A nonvanishing value implies that 
the localized spins contribute to the Fermi surface volume. The second order parameter describes 
singlet-pairing between the localized spins and competes with the Kondo-pairing order parameter. 
Reduction of the carrier density in the conduction band reduces the energy gain due to the formation 
of the large Fermi surface and induces a phase transition to a state with strong singlet correlations 
between the localized spins and a Fermi surface which comprises only the conduction electrons. The 
model thus shows a doping-driven change of its Fermi surface volume. At intermediate doping and 
low temperature there is a phase where both order parameters coexist, which has a gapped large 
Fermi surface and d^2_y2 superconductivity. The theory thus qualitatively reproduces the phase 
diagram of cuprate superconductors. In the second part of the paper we show how the two phases 
with different Fermi surface volume emerge in a strong coupling theory applicable in limit of large 
Kondo exchange. The large-Fermi-surface phase corresponds to a 'vacuum' of localized Kondo sin- 
glets with uniform phase and the quasiparticles are spin- 1/2 charge fiuctuations around this fully 
paired state. In the small-Fermi-surface phase the quasiparticles correspond to propagating Kondo- 
singlets or triplets whereby the phase of a given Kondo-singlet corresponds to its momentum. In 
this picture a phase transition occurs for low filling of the conduction band as well. 

PACS numbers: 71.10.Fd,74.72.-h,71.10.Ay 



I. INTRODUCTION 

The existence and shape of the Fermi surface 
and its change with the hole concentration S in 
the Cu02 planes appears be one of the central is- 
sues in the physics of cuprate superconductors. In 
the overdoped compound Tl2Ba2Cu06+2; magnetore- 
sistance measurements^, angle- resolved photoemis- 
sion spectroscopy (ARPES) 2] and quantum oscillation 
experiments Q show a situation which is reminiscent of 
that in heavy Fermion compounds: despite the participa- 
tion of the strongly correlated Cu3d orbitals in the states 
near the Fermi energy the Fermi surface agrees well with 
LDA band structure calculations which take the Cu3d- 
electrons as itinerant, quantum oscillation experiments 
show the validity of the Fermi liquid description with an 
enhanced band mass. The only moderate mass enhance- 
ment in the cuprates thereby seems natural given the 
large Cu3d-02p exchange constant, W » 1 eV, which 
would give a very high nominal Kondo temperature. 
In the underdoped compounds the situation is more in- 
volved. ARPES shows 'Fermi arcs' Q which however are 
probably just the 'visible' part of hole pockets centered 
near {^,^). This is plausible because the sharp drop 
of the ARPES weight of the quasiparticle band upon 
crossing the noninteracting Fermi surface which must 
be invoked to reconcile the 'Fermi arcs' with the hole 
pocket scenario, is actually well established in insulat- 
ing cuprates such as Sr2Cu202Cl2[5| and Ca2Cu02Cl2[B| 
where this phenomenon has been termed the 'remnant 



Fermi surface'. Meng et al. reported the observation of 
the previously unresolved 'dark side' of the hole pockets 
in underdoped Bi2(Sr2-:rLa:r)Cu06 by ARPES0. Their 
conclusions subsequently were criticized Isj and the issue 
still seems controversial [9j- 

Moreover both the Drude weight in La9_.,.Sr^Cu04[Tol. 
[Tl| and YBa2Cu30y[ll'| as well as the inverse low tem- 
perature Hall constant in La?- r Sr r Cud [lll - [T^ and 
YBa2Cu30j,|ll| scale with S and the inferred band mass 
is constant troughout the underdoped regime and in fact 
even the antiferromagnetic phase lli]. This is exactly 
the behaviour expected for hole pockets. On the other 
hand, for 6 > 0.15 the Hall constant in La2-2:Sr2:Cu04 
changes rapidly, which suggests a change from hole pock- 
ets to a large Fermi surfacefT?]. Quantum-oscillation 
experiments on underdoped YBa9Cu-iOf; sfisl-fisj and 
YBa9Cu/iO«[l9l. [20j show that the Fermi surface has 
a cross section that is comparable to 6/2 rather than 
(1 — (5)/2 as in the ovderdoped compounds. Here it should 
be noted that the mere validity of the Fermi liquid de- 
scription as demonstrated by the quantum oscillations 
is conclusive evidence against the notion of 'Fermi arcs': 
the defining property of a Fermi liquid is the one-to-one 
correspondence of its low-lying states to those of a fic- 
titious system of weakly interacting Fermionic quasipar- 
ticles and the Fermi surface of these quasiparticles is a 
constant energy contour of their dispersion and therefore 
necessarily a closed curve in k-space. On the other hand 
the quantum oscillations cannot be viewed as evidence for 
hole pockets either in that both the Hall constant [21)] and 
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thermopowerf2^ have a sign that would indicate electron 
pockets. Thereby both, the Hall constant and the ther- 
mopower, show a strong temperature dependence and in 
fact a sign change as a function of temperature. This 
sign change is observed only at where superconductiv- 
ity is suppressed by a high magnetic field. At the same 
time neutron scattering experiments on YBa2Cu306.6 in 
the superconducting state show pronounced anisotropy in 
the spin excitations spectrum below 30 meV and at low 
temperatures psj . This indicates an as yet not fully un- 
derstood anisotropic state, possibly to a 'nematic' state 
with inequivalent x- and y-direction in the Cu02 plane. 
Such a nematicity has also been observed in scanning tun- 
neling microscopy experiments on Bi9Sr9CaCu9 0s+^[23| 
and will modify the Fermi surface in some way which may 
explain the unexpected sign. More recently, Sebastian et 
al. concluded from an analysis of the 2"'* harmonic in 
quantum oscillations in underdoped YBa2Cu306+a; that 



the Fermi surface consists only of a single pocket [25 1. 
Since this would rule out the possibility of coexisting 
hole-like and electron-like Fermi surface sheets [1^ and 
since it is hard to imagine that the sole Fermi surface 
sheet of a hole-doped compound is electron-like this re- 
sult would imply that the Fermi surface actually is a hole 
pocket and that the sign of the Hall constant and ther- 
mopower does not reflect the nature of the carriers but 
is determined by some other mechanism. Adopting this 
point of view, the picture of hole pockets centered near 
' ■§) '^ith an area cx 5 would give a simple and consis- 
tent description of ARPES, normal state Drude weight 
and Hall constant and quantum oscillation experiments 
in the underdoped state. Combined with the results for 
the overdoped compounds this would imply that as a 
function of 5 the cuprates undergo a phase transition 
between two states with different Fermi surface volume, 
whereby the Cu 3d electrons contribute to the Fermi sur- 
face volume in the overdoped compounds but 'drop out' 
of the Fermi surface volume in the underdoped regime. 
Interestingly, the superconducting transition in the 
cuprates itself seems to be accompanied by a Fermi sur- 
face change as well. Namely ARPES shows that the 
quasiparticle peaks near (tt, 0), which are hardly distin- 
guishable in the normal state, become very intense and 
sharp in the superconducting state [13, [111. This looks as 
if coherent quasiparticles around (tt, 0) exist only in the 
superconducting state and in fact seem to jurnp into exis- 
tence right at the superconducting transition [23|. A pos- 
sible interpretation would be that the superconducting 
transition occurs between a hole pocket-like Fermi sur- 
face - which does not extend towards (tt, 0) - to a gapped 
large Fermi surface which naturally has some portions 
near (tt, 0). 

Transitions where the correlated electron subsystem con- 
tributes to the Fermi surface volume or not are not unfa- 
miliar in Heavy Fermion compounds. An example is the 
metamagnetic transition in CeRu2Si2 where the Ce 4f- 
electrons, which contribute to the Fermi surface in zero 
magnetic field, seem to drop out of the Fermi volume as 



the magnetic field is increased [3(|. In CeRhi_a; Covins 
the Ce 4f electrons change from localized for x < 0.40 
to itinerant for x > 0.50 as the lattice constant de- 
creases due to substitution of Rh by the smaller Co^Sll. 
Similarly, the localized Ce 4f electrons in CeRh2Si2 _32], 
CeRhlnr, [33j and CelnafsJ] can be made itinerant by 
pressure. It seems that in the case of the metamagnetic 
transition the magnetic field breaks the Kondo singlets 
between Ce 4f spins and conduction electrons whereas in 
the other cases the decrease of the hybridization strength 
between 4f and conduction electrons makes the formation 
of Kondo singlets unfavourable. 

It has been pointed out long ago by Doniach that there 
may be a competition between the Kondo effect and the 
RKKY interaction which then leads to a phase transition 
to a magnetically ordered phase if certain parameters in 
the system are variedfs^. More recently Senthil et al. 
have investigated this question in the context of Heavy 
Fermion compounds and discussed a transition to a mag- 
netically ordered state which is accompanied by a Fermi 
surface transition from large to small as the strength of 
the Kondo coupling is varied 36j. 

In the cuprates one might expect another reason for a 
Fermi surface transition, namely the depletion of the mo- 
bile carriers i.e. holes in 02p orbitals. It is self-evident 
that the gain in energy due to formation of a common 
Fermi sea of mobile 02p holes and localized Cu3d spins 
must tend to zero when the density of mobile carriers van- 
ishes. More precisely, one might expect the band filling 
to play a substantial role when the width of the occu- 
pied part of the conduction band becomes smaller than 
the width of the Kondo resonance at the Fermi level. 
Since the Cu3d-02p exchange constant is large and the 
band filling small this situation may well be realized in 
the cuprates. As the density of holes in 02p orbitals is 
reduced one would thus expect that at some point the 
localized spins drop out of the Fermi surface so as to op- 
timize their mutual superexchange energy instead. Due 
to the near-two-dimensionality of the cuprates, however, 
the transition is not to a magnetically ordered state, but 
to a 'spin liquid' with strong nearest-neighbor singlet cor- 
relations instead of true antiferromagnetic order. 
Cuprate superconductors are frequently described by a 
single band Hubbard model or the t-J model. Whereas a 
Kondo-lattice-like Hamiltonian can be derived by lowest 
order canonical perturbation model from the so-called d- 
p model for the Cu02 plane [s^ these single band models 
are obtained in a subsequent step which is valid in the 
limit of large Kondo coupling between 02p and Cu3d 
electrons, so that the Kondo singlet (Zhang- Rice- singlet) 
extends over little more than only one plaquette. Yet, 
these models should show a Fermi surface transition as 
well if they are equivalent to the Cu02 plane. Exact di- 
agonalization studies of the t-J model have indeed shown 
that the Fermi surface at hole dopings < 15% takes 
the form of hole pockets (38l - l40j . that the quasiparticles 
have the character of strongly renormalized spin polarons 
throughout this doping range (4l| - |43j and that the low en- 
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ergy spectrum at these doping levels can be described as 
a Fermi liquid of spin 1/2 quasiparticles corresponding to 
the doped holes [4J|. A comparison of the dynamical spin 
and density correlation function at low[45il46| {S < 15%) 
and intermediate and high = 30 — 50%) hole dop- 
ing moreover indicates [47] that around optimal doping 
a phase transition to a state with large Fermi surface 
takes place in the t-J model. A study of the electronic 
self-energy in the single-band Hubbard model has indi- 
cated that there such a transition takes place as well[48|. 
Contrary to widespread belief such hole pockets can be 
completely consistent with the Luttinger theorem for the 
single band Hubbard model[48l]. 

To study the issue of the Fermi surface transition fur- 
ther we performed a Hartree-Fock treatment of a two- 
dimensional Kondo lattice with an additional Heisenberg 
exchaneg between localized spins. This is presented in 
Section H. Interestingly it turns out that such a transition 
not only exists but is generically accompanied by (1^2 _y2 
superconductivity in the mobile carrier system which 
might provide an explanation for the phase diagram of 
the cuprates. In simplest terms superconductivity occurs 
because the coherent Kondo-pairing between d-spins and 
c-electrons 'transfers' the singlet pairing between the lo- 
calized d-spins to the mobile c-electron sytem. 
Since a mean-field treatment may not really be expected 
to be valid in the limit of large superexchange W between 
localized and conduction electrons, which is the region of 
physical interest for the cuprates, we also show how the 
two phases with small and large Fermi surface may be 
understood in the limit of large W. The small Fermi 
surface phase is discussed in section HI and is similar 
to the lightly doped Mott insulator as discussed in Ref. 
[49} . The phase with large Fermi surface is discussed in 
section IV, closely following Refs. [50,] and [51,]. Section 
V discusses the possibility of a phase transition between 
small and large fermi surface within the strong coupling 
theory. Section VI gives a summary and discussion. 



II. MEAN FIELD THEORY 

We study a Kondo lattice model which consists of a sin- 
gle metallic conduction band - described by the Fermionic 
operators c- ^ - and a lattice of localized spins, described 

by dl^. The lattice sites - labelled by i - form a simple 
cubic lattice and there is one localized spin and one s- 
like conduction orbital in each unit cell. A similar model 
with two 02p-like conduction orbitals per unit cell can 
be derived by canonical transformation for the Cu02 
plane^J. We augment the model by a direct Heisen- 
berg exchange between nearest neighbor localized spins. 
A similar model for the Cu02 plane has been studied 
previously [52|. To make more easy contact with the 
CUO2 planes we consider the operators cl ^ and dl ^ to 
create holes rather than electrons. We denote the density 
of holes in the conduction band hy 6 - the total densisty 



of holes/unit cell thus is 1 -(- 5. The Hamiltonian reads 
H = Ht + Hw + Hj, 



H.I 



(1) 



where 



2 '^i,a'^a,l3^i,l3i 



(2) 

with a the vector of the Pauli matrices and analogous 
definitions hold for Sd.i and n^.i- The model is to be 
considered in the sector of the Hilbert space where all 
d-orbitals are singly occupied. 

In a first step we drop the d — d exchange Hj. We use 
the identity 



Sd,' 



s] = 



d^ 



d^ 



and apply the Hartee-Fock approximation 



(3) 



(4) 



Whereas the orginal Hamiltonian conserves the number 
of d-holes at each site, this does not hold true for the 
mean-field Hamiltonian - which is clearly a severe draw- 
back of the theory. Although the present decoupling 
scheme is different from that used by Senthil et al. [36l] 
the theories can be converted into each in the case J = 
by performing a particle-hole transformation for the d- 
electrons. 

There are two side conditions to be obeyed: one for the 
total hole number, Nh, the other one for the number of 
d-holes, Nd, which must be equal to N, the number of 
d-sites in the system. We enforce these by Lagrange mul- 
tipliers i-j, and A, respectively, and, introducing the vec- 



tor <.<T - v-k,^' 
Hamiltonian 



cj^ o-j c?k.ff), we obtain the Fourier transformed 



H-nNh- XNd 



E 

k,CT 



""La Hk,a-"k,a 



-2Ned + 2N^i + N 



2W' 



(5) 

where Nh [Nd) are the operators for the total number of 
holes (number of d-like holes), N is the number of sites, 
and 



Hk,cr 



ek - M 



-sign(a)- 



X + 2/1. 



(6) 
(7) 
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The parameter Acrf describes coherent singlet formation 
between the conduction electrons and the localized spins. 
The essence of the Kondo-effect is the quenching of the 
magnetic moments that means a localized spin at site i 
forms a singlet with a conduction electron (this would 
be the Zhang-Rice singlet in the case of cuprate super- 
conductors). The phase of this local singlet then is in 
principle undetermined and the above mean-field decou- 
pling describes a state, where this phase is uniform over 
the whole system. 

Next, being the expectation value of two cre- 
ation/annihilation operators. Ad has some similarity to 
a superconducting order parameter. Since the pairing is 
not between time-reversed states, however, the resulting 
gound state is not superconducting. This will be appar- 
ent from the fact that it has a well-defined Fermi surface. 
The Hamiltonian can be diagonalized by the transforma- 
tion 



7k.2,cr 



sign{a)vi^d^i,^, 

and we obtain the two quasiparticle bands 
1 



(8) 



(ek + ± Wi^) - fj,, 



Ek - Sd 

2VFk 

Ek - 

2Wu 



1/2 



1/2 



(9) 



In the equation for Uk it is assumed that A^^ > 0. The 
self-consistency equation for Acd becomes 



1 



W 

TV 



E- 



sinh 



(^) 



k ^kcosh(^) +cosh(/3( 
Next, it is straightforward to show that 



2 



(10) 



27V-iVrf + 7V, = ^^7i_,^^7k,,,, (11) 



k.cr u—1 



where Nd {Nc) are the operators for the number of c-like 
(d-like) holes. If we demand that Nd — N the operator 
Nh of the total hole number becomes 



(12) 



k,cr l'=l 



so that we have a 'large' Fermi surface which comprises 
the localized d-holes. This, however, will hold only if 
there is exactly one localized spin/d-site. When Acd is 
zero the c-holes and d-holes are decoupled and the con- 
dition Nd = N puts a dispersionless cZ-band right at fi, 



0.14 
0.12 
0.1 
0.08 
0.06 
0.04 
0.02 



W/t=2.0 






W/t=1.8 — - 






■ W/t=1.6 - 






W/t=1.4 






- W/t=1.2 


- 




W/t=1.0 - - - 








-I 






FIG. 1: Critical temperature for the onset of the parame- 
ter Acd with the true 2-dimensional nearest-neighbor-hopping 
band structure (top) and a constant density of states in the 
Intervall [-At : At] (bottom). 



i.e. ed = A*. The Fermi surface then is that of the de- 
coupled c-holes and accordingly has a volume which does 
not comprise the d-holes. 

As already stated for Acd —J' we have — > /i and insert- 
ing this into (|10|) we obtain the equation for the critical 
temperature 



W 

2iv4^ek 

k 



1 



■ tanh 



(13) 



Figure [T] shows Tc as a function of S. This is shown 
both for the true 2-dimensional nearest-neighbor-hopping 
dispersion and for a conduction band with a constant 
density of states in the interval [—At : At]. It is obvious 
that Tc —> as 5 —>■ 0. In the limit of small Tc we can 
obtain a rough approximation by replacing the integrand 
in (|T3| by (|ek — + Tc/2)~^. For the constant density 
of states we find in this way 



T, = 2t e-8*/w'^5(2-5). 



(14) 



As expected, the energy gain due to the formation of 
common large Fermi surface thus goes to zero when the 
density of mobile carriers vanishes. For sufficiently low 
temperature the large Fermi surface thus is formed for 
any carrier concentration, but if there is a competing 
term in the Hamiltonian there may be a phase transition 
at finite doping. 

We now introduce the d-d Heisenberg exchange cx J. 
We decouple the d-d Heisenberg exchange in the same 
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way as the c-d exchange: 

nd.incu 1 



Sd.i ■ Sd 



J - J 



dl^d] , -d] ,d 



t 



Inclusion of this term doubles the dimension of 
the matrices to be considered. Introducing Vk = 
(cjj. ^, d jj. ^, c_jj. ^, djj. ^) the Hamilton matrix becomes 



V 



2 








-ek + M 

A!. 



-ArfdTk 

Acd 
2 

-<^d + ) 



1 



7k = - (cos(fc^) ± cos(fcj/)) 
Add 



2 



(15) 



The two parameters, Acd and A^d now have to be deter- 
mined self-consistently. Unlike the on-site order param- 
eter Acd, the d — d pairing amplitude is a 'bond related' 
quantity and thus may have different sign for bonds along 
X and y so that we may have s-like pairing or d-like pair- 
ing. These two possibilities have to be considered sepa- 
rately. 

Before studying the full problem we briefly consider the 
case Acd = 0. In this case the d- and c-holes are again de- 
coupled and we only have to treat the d-electron system. 
The chemical potential is determined by the conduc- 
tion holes alone and we must have td = [i- The result 
will not depend on 5 or 14^ and the temperature depen- 
dence is universal if temperature is measured in units of 
J. Finally, there is no difference between 7s_d(k) because 
cos(fcj,) = — cos(fcj, -f tt) so that the d-like pairing simply 
corresponds to a shift of the Brillouin zone by (tt, 0). The 
self-consistency equation reads 



2J 



sinh(/37(k)Add) 
+ cosh(/37(k)Add) 



(16) 



and the temperature for the phase transition is — J /A. 
The gap at T = is found to be 



Add = jY.\^m 

k 

= J -0.81. 



(17) 



The localized spins aquire a nonvanishing dispersion, 
which is unphysical and an artefact of the mean-field ap- 
proximation. The problem is somewhat lessened in that 
the localized electrons at least have no Fermi surface. To 
see this we note first that the dispersion consists of two 
bands Cd ± Add7(k). Since ed = the lower of these is 
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FIG. 2: Development of the two order parameters with in- 
creasing J. The other values are W/t — 1.6, 5 = 0.2 and 
T/t = 0.001. 



completely filled. The momentum distribution function 
for the d-electrons then becomes 



nd.k 



1 



(18) 



so that at least the momentum distribution of the d- 
electrons is consistent with localized electrons. 
Next we switch to the full problem of two coupled order 
parameters i.e. Acd ^ and Add 7^ 0. All results pre- 
sented below have been obtained by numerical solution 
of the self-consistency equations on a 400 lattice with pe- 
riodic boundary conditions. Study of the variation with 
lattice size shows that this implies a reasonable conver- 
gence. 



As already mentioned Add may be s-like and d^: 



-like. 



This difference will now matter because the relative po- 
sition of the lines of zeroes in the rf-electron dispersion 
and the Fermi surface of the c-electrons makes a physi- 
cal difference. To decide which of the two symmetries is 
realized we consider the free energy/site given by 



/3 



k 1^=1 



A? A^ 

cd ^ dd , / £ 



1)a^- (19) 



Numerical evaluation shows that the s-like state is never 
realized. 

Figure [2] then shows an example of the development of 
the two coupled order parameters as J is switched on. It 
is quite obvious that increasing Add is detrimental to Acd 
i. e. the two order parameters are competing. Nonvan- 
ishing J thus may introduce a phase transition at T = 
as a function of doping between the two phases. 
To discuss the phase diagram we note first that we can 
distinguish two regimes: at low doping the critical tem- 
perature Tc for the onset of Acd will be below J/4, which 
is the critical temperature for the onset of Add. This 
means that as the temperature is lowered, a nonvanish- 
ing Add sets in first and Acd then must develop on the 
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FIG. 3: Temperature dependence of the two order parameters 
for Wjt = 1.6, J It = 0.05, (5 = 0.32. 
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FIG. 4: Self-consistent solution for function of 

temperature for various dopings 5. The other values are 
Wjt = 1.6, J It = 0.05. 



'background' of this nonvanishing I^dd- When Tc > J/4 
on the other hand, we first have a nonvanishing Acd and 
Atjrf develops at lower temperature. As will be seen next, 
the two regimes are more different than might be ex- 
pected at first sight. For the time being we fix W jt = 1.6 
and Jjt = 0.05. As can be seen in Figure [H the doping 
where = J/4 then is approximately i5 = 0.3. Figure [3] 
then shows the self-consistent A's as a function of tem- 
perature for b — 0.32. Both parameters show a fairly 
conventional behaviour at the two phase transitions with 
the characteristic \JTc — T behaviour and the same is 
seen whenever A^d sets in at higher temperature. 
The situation is quite different for 5 < 0.3 as can be 
seen in Figure HI For most dopings there is now a 
small but finite temperature range where two solutions 
(Acd ^ 0, Add 7^ 0) exist (in addition there are the so- 
lutions with Add = and Acd — 0). Calculation of the 
free energy shows, that it is always the solution with the 
higher Acd which has the lower free energy, i.e. this is 
the physical state. To clarify the nature of the phase 
transition we fix Acd, determine all other parameters 
(/i, Cd, Add) self-consistently and evaluate the free energy 
/ as a function of Acd- Points on the resulting curve 
/(Acd) which are stationary with respect to variations of 
Acd are stationary with respect to variations of all param- 
eters and therefore are solutions to the self-consistency 
equations. The result is shown in Figure [S] At high 
temperature /(Acd) has only one extremum, namely a 
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FIG. 5: Top: Self-consistent solution for Acd as a function of 
temperature for 5 — 0.25, see Figure |4l Bottom: Scans of the 
free energy as a function of Ad whereby all other mean-field 
parameters have been obtained self-consistently. 



minimum at Acd = 0. As the temperature is lowered, 
however, one can recognize a 'wiggle' in the curve which 
develops into a maximum and a minimum. These corre- 
spond to the two solutions with nonvanishing Acd. Next 
there occurs a level crossing between the two minima - we 
thus have a first order transition - and as the temperature 
is lowered further the maximum for Acd 'absorbs' the 
minimum at Acd = 0. From then on we have only the 
minimum with Acd 7^ and the maximum at Acd = 0. 
This behaviour can be seen at almost all dopings below 
S. It is only very close to 6 that there is only one solu- 
tion but the Acd(r) curve is very steep as Acd 0. The 
important finding then is that the temperature-induced 
transition between small and large Fermi surface is a 
order transition in this doping range. This is not surpris- 
ing in that a 2"'' order transition involving competing 
order parameters may become 1"* order js^. 

Next we consider the dependence on hole doping 6, 
in particular the question if the Fermi surface transition 
now occurs at finite doping even at zero temperature. 
Figure [H which shows the self consistent solutions for 
Acd at T = for different W/t and J/t demonstrates 
that the answer depends on the magnitude of J/t. For 
small J a nonvanishing solution exists down to i5 = 0. In 
fact there are two such solutions, but computation of the 
free energy shows that only the solution with the larger 
Acd is a minimum and moreover has a lower free energy 
than the solution with Acd = 0. Accordingly, there is 
no phase transition down to (5 = for small J/t. If J 
increases, however, there is now an extended range of 
small i5 where no solution with Acd 7^ exists. As can 
be seen in the Figure [H] this change occurs in that the 
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FIG. 6: Self-consistent solutions for Acd as a function of dop- 
ing for T = and different J. The value Wjt = 2 {Wjt = 1.6) 
in the upper (lower) part of the figure. 



maximum and minimum of /(Ac^) 'merge' for some 8. 
The hourglass shape formed by the two Ac(i((5) curves 
immediately below the critical J jt can be clearly seen in 
the upper part of the Figure. It has been verified, how- 
ever, that for large values of J ji where the transition has 
occured there is no more solution for low 5. This means 
that for larger J/t a phase transition does occur and there 
is a considerable range of i5 where the mean-field theory 
predicts a spin-liquid (i.e. A^^ ^ 0) with a small Fermi 
surface. 

We proceed to a discussion of the phase diagram whereby 
we consider the more interesting case of larger J jt. An 
example is shown in Figure [71 the band structures for the 
various phases are shown in Figures [S] and [HI To begin 
with, we can distinguish four phases. At high temper- 
ature (phase 1) both self-consistent parameters are zero 
which implies that the conduction holes are decoupled 
from the localized spins and the localized spins them- 
selves are uncorrelated. At higher doping and low tem- 
perature there is a phase with Acd ^ 0, Acd = (phase 
2). The band structure (see Figure [51 Top part) shows 
that there is a large Fermi surface with an enhanced band 
mass. This phase most likely corresponds to the over- 
doped regime in cuprate superconductors. 
For low doping and low temperature, on the other hand, 
there is a phase with A^j^ 7^ 0, Acrf (phase 3). 
This has strong singlet correlations between the local- 
ized spins and a small Fermi surface which is centered at 
(tTj-tt) in mean- field theory, (see Figure [H Bottom part). 
This phase likely corresponds to the 'pseudogap' phase 
in cuprate superconductors. Thereby the strongest defi- 
ciency of the mean-field theory consists in the neglect of 
any correlation between the localized spins and the con- 
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FIG. 7: Top: Mean-field phase diagram for the parameter 
values Wjt = 2, J jt = 0.325. 

Bottom: Difference in entropy/site between the phases 3 and 
4. 




FIG. 8: Top: Band structure for phase 2 (T = 0.005, 5 = 

0. 60). The two vertical lines at /i give the Fermi momenta, 

1. e. the crossings of the quasiparticle band through [s,. 
Bottom: Band structure for phase 3 (T = 0.005, A = 0.25). 
Other parameter values are W jt = 2, J /t = 0.325. 



duction holes. One may expect that coupling of the con- 
duction hole pocket around (tt, tt) to the antiferromag- 
netic spin fluctuations of the localized spins will create 
hole pockets centered near (^, ^). In any way, however, 
the Fermi surface does not comprise the d-electrons in 
this phase. 

Finally, at intermediate doping and low temperature 
there is a 'dome'-like region in the phase diagram in 
which both self-consistent parameters are different from 
zero (phase 4). Inspection of the band structure in Figure 
[51 shows that phase 4 does not have a Fermi surface, but 
rather a gap with a node along the (1, 1) direction, i.e. 
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FIG. 10: Pairing amplitude for conduction holes Acc as a 
function of temperature, W jt = 2.0, J /t = 0.325. 



FIG. 9: Band structure for phase 4 (T = 0.005, 5 = 0.45). 
Entire bandwidth (top) and closeup of the region near (bot- 
tom). Other parameter values are W/t — 2, J/t = 0.325. 

the band structure expected for a superconductor with 
da;2_j,2 order parameter. The Fermi momentum along 
the (1,1) direction coincides with that of the large Fermi 
surface, however, so that we have a gapped large Fermi 
surface in this phase. Across the phase transition on 
the low doping part of the dome the Fermi surface thus 
changes from a hole pocket centered on (tt, tt) to a gapped 
large Fermi surface, whereas the transition on the high 
doping side of the dome corresponds to the opening of 
a d^2_j,2-like gap on the large Fermi surface. Accord- 
ingly, we consider the question whether phase 4 is super- 
conducting. To that end we consider the c-like pairing 
amplitude 

^ ^ E^(l^)(4T^-ki + C-H^kt)- (20) 

k 

Figure [10] shows that this is indeed different from zero 
within the phase 4. Interestingly the Hamiltonian does 
not contain any attractive interaction between the con- 
duction holes. Rather, the singlet pairing between the lo- 
calized spins as described by the parameter Add is 'trans- 
ferred' to the mobile conduction hole system by the co- 
herent Kondo pairing amplitude A.cd- 
Figure [11] shows the gap and the c-like spectral weight for 
several momenta along the (1, 0)-direction and demon- 
strates that the present theory qualitatively repro- 
duces a well-known anomaly in cuprate superconductors: 
whereas the gap size near (tt, 0) increases with decreasing 
Tc - in contrast to what one expects from BCS theory - 
the spectral weight decreases to zero 'HI]. At least the 
behaviour in Fig. [TT]is easy to understand: the gap size 
around (tt, 0) is determined by Add, the c-like spectral 
weight by Acd which governs the degree of mixing be- 
tween c- and d-like bands. Lower Tc implies a lower value 
of Acd at low T and - since the two parameters compete 
with each other - a larger value of Add- To conclude this 
section we briefly summarize the results of the the mean- 
field theory: The competition between Kondo-coupling 
and Heisenberg-exchange between localized spins leads 
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FIG. 11: Energy gap (top) and spectral weight of the band 
below fi (bottom) as a function of doping for some momenta 
along (1,0). Parameter values are W/t = 2.0, J/t = 0.325, 
T/t = 0.05. 

to a doping-driven phase transition between states with 
different Fermi surface volume. The phase for high dop- 
ing is characterized by a complex order parameter which 
describes coherent singlet-pairing between localized and 
conduction holes and leads to a Fermi surface volume 
that includes the localized spins. The low doping phase is 
characterized by an order parameter which describes sin- 
glet correlations between the localized spins and there is 
no coupling between conduction holes and localized spins. 
The Fermi surface thus is a hole pocket with a volume 
which does not include the localized spins. Interestingly 
there is an intermediate phase where both order param- 
eters coexist. This phase has a gapped large Fermi sur- 
face and shows superconducting correlations between the 
conduction holes. These may be interpreted as the sin- 
glet correlations between the localized spins being 'trans- 
ferred' to the conduction electrons by the coherent Kondo 
singlet formation. This implies that the superconducting 
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transition has a very different cfiaracter on the under- 
doped and overdoped side of the superconducting dome: 
whereas on the overdoped side the transition looks fairly 
conventional, with a gap operning in the large Fermi sur- 
face, on the underdoped side the transition corresponds 
to the sudden emergence of a gapped large Fermi surface 
from the high-temperature phase with a hole pocket. A 
somewhat discomforting feature of the transition in the 
underdoped range is the fact that it is 1^* order. Lastly 
we mention that the criticial temperatures and dopings in 
the present study are very different from the cuprates but 
this can be hardly a surpise because we are studying a 
different model (single band Kondo lattice rather than 
Kondo-Heisenberg model) and mean-field theories can 
not be expcted to produce quantitatively correct tran- 
sition temperatures anyway. 

III. STRONG COUPLING THEORY - SMALL 
FERMI SURFACE 

The mean-field theory in the preceding section may be 
expected to be good for t:^ W,J although even there re- 
laxing the constraint on the occupation of the d-orbital is 
problematic. The phsical regime of parameters, however, 
is rather W ^ t J. In fact a considerable deficiency 
of the mean-field treatment lies the fact that in the hole 
pocket phase the d- and c-holes are completely decou- 
pled from one another, which clearly is unrealistic for 
large W/t. In the following sections we therefore give a 
description of the two phases in a strong-coupling picture 
which may be expected to hold best in the limit W/ 1 ^ 1 . 
The discussion of the small-Fermi-surface phase in the 
present section thereby closely follows the theory for the 
lightly doped Mott insulator given in Ref. [ioj . 
The essence of the strong coupling theory is the approx- 
imate diagonalization of the Hamiltonian in a suitably 
chosen truncated Hilbert space. To construct this trun- 
cated Hilbert space for the small-Fermi surface state we 
start from the case (5 = and consider a state |^o) of 
the d-spin system which has exactly one electron/site, is 
invariant under point group operations, has momentum 
zero and is a spin singlet. These are the quantum num- 
bers of a vacuum state and indeed jvPo) will play the role 
of the vacuum state in our analysis. The only property 
of l^'o) which is relevant for the quasiparticle dispersion 
and total energy is the static spin correlation function 

X., -(*o|S.-S,|vI/o). (21) 

We consider Xij as a given input parameter. We assume 
it to be antiferromagnetic and of short range i.e. 

X., =Coe''2-(R.-R.)e-^^^ (22) 

where Q = (tt, tt). A more detailed discussion is given in 
Ref. It will be seen below that within the frame- 

work of the present theory only the nearest neighbor spin 
correlation xio has any relevance for the results. 



Next, we define the following operators which add a con- 
duction hole to the system: 

4t - 7|(4t^4~4,;S+), 

^U.T = 4Av (23) 

Here capital/small letters are used for local- 
ized/conduction holes, Ci^a = Ci,CT(l — rig), and 
Ni^a = -/Vv(t(1 — Na). The operators ([23|) add a conduc- 
tion hole to the system, thereby change the z-component 
of the total spin by +1/2 and produce either a local 
singlet (al^) or one of the two components of a local 
triplet (bli^ and bl2^)- Analogous operators which 
change the z-componcnt of the total spin by — ^ are 
easily constructed. 

Using these operators we can now write down the basis 
states of the truncated Hilbert as 




(24) 

Thereby it is understood that all sites (v, j^, k\) are pair- 
wise different from each other, that means no two oper- 
ators are allowed to act on the same site. The states 
([24l) have singlets and triplets at specified positions and 
we will treat these as spin-^ Fermions with a hard-core 
constraint, which is the key approximation of the theory. 
Fermions are the only meaningful description for these 
particles because operators of the type which refer 

to different sites anticommute. Since (^o|ai,t'^I I'^'o) = 

(^'o|^»4,T^Ii,tl*o) = (^'o|6z,2,t^I,2,tl*o) = 5 the states 
(|24p are approximately normalized. The issue of the nor- 
malization of the states has been discussed in detail in 
Ref. ^1 where it was concluded that this normalization 
will be a good approximation in the limit of short spin 
correlation length that means the case of a 'spin liq- 
uid' which is the case of interest e.g. in the underdoped 
cuprates. A more detailed discussion of this issue and 
others is given in Appendix A. 

As already stated we treat the singlets and triplets as 
Fermionic quasiparticles i.e. the states (1^^ are repre- 
sented by states of Fermionic spin-^ quasiparticles 

(n«i..] (n/i-.) (n/L^^^. ) i^o). 

(25) 

Operators in the quasiparticle Hilbert space then are de- 
fined by demanding that their matrix elements between 
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the states ((25)) are identical to those of the physical op- 
erators between the corresponding states ([M)) . 
We illustrate this by setting up the quasiparticle Hamil- 
tonian. Straightforward calculation gives the following 
matrix elements in the physical Hilbert space: 



the physical Hilbert space: 



{■^o\bj,i.,t Ht 4tl*o> 



1 , X±j 
2 

i _ 
8 6 

1 X 



1 , Xy 



{^oK2.^ Ht &I,2,tl*o) 



3 ' 



(26) 



They represent the decomposition of hopping into events 
of annihilation and creation of singlets and triplets. The 
hopping term Ht for the quasiaparticles thus takes the 
form 



(27) 



where 



and 



Tij — tij 



^la = («L'^ll,<T'^I,2,a)> 



(28) 



4 ^ 4 4 "I" ) 

V2(i - ^) 72(1 + 2^) 



(29) 

Thereby we have to keep in mind the hard-core constraint 
between the quasiparticles. In addition it has been as- 
sumed that 'nearby' quasiparticles do not modify the ma- 
trix elements describing the propagation of a given quasi- 
particle substantially. Again, this will be justified in the 
limit of short spin correlation length C,. 
The key simplification is that in our restricted Hilbert 
space the Kondo-coupling term Hw takes the simple 
form 



3VF + W^V^V^it 7 l^r^\ 

Hw = — "''^"^'^ + T Kn,A,n,a- (30) 



In the strong coupling limit W/t 1 we thus expect to 
obtain something like one lower and two upper 'Hubbard 
bands' separated by an energy of order W . 
It remains to represent Hj in the new basis. Straightfor- 
ward computation gives the following matrix elements in 



(5'nld,- Si 

\ ^ U 1 '^l. V ^4 


• S,) aU 


1 ^ u/ 


= 


(*ol&.,i.T (S. 




l*o) 


_ Xij 

~ 6 


(*o|&^,2.T (S. 


• s,) 4^1 


l*o) 


_ Xij 






l*o) 


= 


(*o|fo^,2,t (S» • 


s.) 


l*o) 


3^2 






l*o) 


Xii 
6 ■ 



so that the d — d exchange for the quasiparticles takes 
the form 



Hj = zJxw E E vJ_<^Xv.,,, 



zNJxi 



(32) 



where xio denotes the nearest neighbor spin correlation 
function and 



K 



_1 1 VI 

-^3 3 

1 _1 \/l 

3 3 

\/2 \/2 _ 2 

3 3 3 



(33) 



The additional constant is the constribution from the 
'spin background' whereby the correction due to the z 
broken bonds/quasiparticle is accounted for in the first 
term. 

The exchange term also has matrix elements of the type 

{^o\ki,Ai,t (S. • S,) a^Xrl^o) = ^ + if ■ (34) 

Their contribution to the total energy will be cx 6^ and 
we will neglect these. 

Next we consider the hole count. The number of local- 
ized holes always is N , the number of sites in the system, 
whereas a self-evident expression for the number of con- 
duction holes is 

= E («L«k,. + &Li,.^k,i,. + ^L2..^k,2..) • (35) 

k,CT 

This will give rise to a Fermi surface whose volume cor- 
responds to the doped holes but does not include the 
(i-electrons. 

The resulting quasiparticle band structure is shown in 
Fig. [12] For a low density of quasiparticles it will be a 
reasonable approximation to neglect the hard-core con- 
straint, because the probability that two particles occupy 
the same site and thus violate the constraint is small. The 
Fermi surface then takes the form of a hole pocket cen- 
tered at Q = (7r,7r). For a semi-quantitative discussion 
we may use 2""^ order perturbation theory for the disper- 
sion ei^k of the lowest quasiparticle band around Q. This 
gives 
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FIG. 12: Dispersion of the quasiparticle bands for the phase 
with small Fermi surface. Parameter values are W/t — 4, 
J/t = 4, xio = —0.28. The dashed line is the dispersion of 
the lowest band as obtained by 2"'' order perturbation theory. 



The band minimum is at (0,0) for xiQ > ~j at 
(tt, tt) otherwise. The third term on the r.h.s. favours 
a xio which is positive, the contribution from the 'spin 
background' favours a negative xio- 

Introducing k = k— Q or k = k, depending on whether 
the minimum of the dispersion is at k = (tt, tt) or k = 
(0, 0) we find 



ei.k = 



Co 



Ci = ±t 



Co + ClK 

1 

4 + 





A8t^ 








- W 


a- 


3 / 



Xio 



24i2 



Xw 
3 



and the ground state energy/site becomes 
Co = cqS + cittJ^ + 2 Jxi 
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(37) 



(38) 



This can now be used to search for a minimum of cq as 
a function of xw- It turns out that for most parameter 
values this expression favours a xio which is as negative 
as possible. We thus set xio = —0.33, the value reahzed 
in the ground state of the Heisenberg antiferromagnet. 
This is obviously an absolute lower bound for the nearest 
neighbor spin correlation function that can be realized by 
any wave function. 



IV. STRONG COUPLING THEORY 
FERMI SURFACE 



LARGE 



The strong coupling description of the phase with a 
large Fermi surface has been given in Refs. [13, HH - see 
also Ref. 56] for a different derivation - and here we sketch 
the derivation only roughly. We start again by defining 
the 'vacuum state' for the theory which now reads 



l^-o) 



-N/2 



ft 



)|0). (39) 



This is a product of local singlets and is the ground state 
of the model for t/W = J/W = and a 'hole doping' 



of (5 = 1. Acting with the hopping term Ht onto p9|) 
produces charge fiuctuations, i.e. states of the type 



(40) 



In this state both cells, i and j, have a total spin of 1/2 
which is carried by the unpaired d-hole-spin. We now 
identify the quasiparticle-states of a single unit cell i as 
follows: 



|0) ^ ^{cld] 



1 

ayo) ^ 4^10) 
«Ulo) ^ 4,10) 



)|0 



J „t jt 



|0) 



&UlO) ^ 4ct,4j0) 
In this representation, the state (HUl) becomes 



sign{a) b]^^ ajjO). 



(41) 



(42) 



Just as ([39|). the state (j40|) is an eigenstate of Hw with 
eigenvalue —{N — 2)-^. To keep track of this large en- 
ergy change we ascribe an energy of iW/A to each of the 
quasiparticles so that the representation of Hw becomes 



V- r t 



(43) 



The particles are created/ annihilated in pairs and by sub- 
sequent application of the hopping term they also can 
propagate individually. The procedure to be followed 
then is entirely analogous as above: we consider a quasi- 
particle Hilbert space whose basis is formed by states of 
the type 



n 




(44) 



As was the case for the small-Fermi-surface phase we as- 
sume that the quasiparticles obey a hard-core constraint, 
i.e. all sites ii, and are pairwise different from each 
other. The corresponding states in the Hilbert space of 
the physical Kondo-lattice are 



2-(N-N^-N,)/2 



Ni 



n 




where 



i/=i 



(46) 



Again we construct operators for the quasiparticles by 
demanding that the matrix elements of operators be- 
tween the quasiparticle states (|44l) are equal to those of 
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the physical operators between the corresponding Kondo- 
lattice states (H5|) . We thus obtain the Hamiltonian (see 
Ref. [51] for details) 

Iv-r/ t , 31^.. ^ . 

_ 1 ^ (6^ali^__ + i/.c). (47) 

where ek = 2i(cos(fc2;) + cos(fcj^)) denotes the dispersion 
relation for the conduction holes. If we assume that the 
density of quasiparticles is low - which will hold true in 
the limit of large W/t and close to (5 = 1 - it will again 
be a reasonable approximation to relax the hard core 
constraint. The Hamiltonian then can be diagonalized 
by the ansatz 



(48) 



and, introducing A=3VF/2, we obtain the quasiparticle 
dispersion 



±A/e?. + A2 



~ 4 2 4A' 



(49) 



where the second line holds in the limit W/t ^ 1. The 
operator of total hole number is given by 



k.cr 



(50) 



The first line follows from the fact that the vacuum 
state (15^ contributes 2N holes, and that each hole- 
like/electron-like quasiparticle increases/decreases the 
number of holes by one. The second line shows that 
the lower of the two quasiparticle bands is filled such 
that the Fermi surface has a volume which corresponds 
to conduction holes and localized spins, i.e. this state 
has a large Fermi surface. The apparent contribution of 
the localized spins to the Fermi surface can be under- 
stood as follows: at a conduction hole density of 1/unit 
cell the state resulting from the above construction has 
a hole number of 2 /unit cell. At the same time the state 
has an energy gap of order 3M^/2, the energy required to 
break two singlets. As far as the hole number and Fermi 
surface - or rather: absence of a Fermi surface - are con- 
cerned, this state therefore is completely equivalent to a 
band insulator, provided the localized spins contribute to 
the total hole number. And since the quasiparticles in- 
troduced by changing 5 are spin-1/2 Fermions the Fermi 
surface at lower conduction hole density is the same as 
that of a 'doped band insulator' and the localized spins 



apparently contribute to the Fermi surface volume. This 
is therefeore simply the consequence of the fact that the 
quasiparticles are spin-1/2 Fermions. 
A notable feature of the above theory is, that it actually 
incorporates the kind of broken gauge symmetry which 
became apparent already in the mean-field treatment: 
all single cell basis states of the phsical system - i.e. the 
states on the r.h.s. of dH]) - are defined only up to a 
phase factor. The quasiparticle Hamiltonian (|47|) then 
holds true only if the phases of all states in a given unit 
cell are equal so that no net phase enters the Hamilto- 
nian. However, one might choose e.g. the phase of the 
states corresponding to &|^^|0) and aj^|0) equal to unity 
and the phase of the singlet state corrsponding to |0) to 
be exp(i(f)j). In this case all matrix elements in the real- 
space version of (|T7)) would aquire extra phase factors. 
Such a phase-disordered state might be adequate to de- 
scribe the Kondo lattice at temperatures lower than the 
Kondo temperature but higher than the coherence tem- 
perature. In the mean-field description we had a conden- 
sate of singlets described by the constant order parameter 
Acd - since the constraint of localized d-holes cannot be 
taken into account rigorously in mean-field theory this is 
probably the best approximation to an array of phase- 
coherent Kondo singlets. 

It remains to discuss the direct d—d exchange oc J. In the 
vacuum state the d—d exchange can promote two singlets 
on neighboring sites into triplets. The matrix element for 
this transition is J/2 so that second order perturbation 
theory gives the correction to the energy/site 



2 8W 



(51) 



This is cx J'^ /W and thus much smaller than the direct 
d—d exchange cx J in the phase with small Fermi surface 
(see equation ((32)) ). This term will be even less important 
because this contribution occurs only if both cells are 
unoccupied by quasiparticles. 

In addition to this contribution of the vacuum there is 
also a contribution of the quasiparticles to the exchange 
energy because the spin of the quasiparticles is carried 
by the d-electron. For example we have 



(52) 



which again can be verified by comparing matrix ele- 
ments of the r.h.s. and l.h.s. between quasiparticle states 
and physical states of the Kondo lattice. 
If we relax the hardcore constraint a straightforward cal- 
culation then gives the contribution to the ground state 
expectation value/site from the exchange between quasi- 
particles of 



(Hj) 
N 



= -3J[ 



2 |^^E^kWkWk/l,kj ] 



(53) 
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where /i^k denotes the ground state occupation number 
of the lower quasiparticle band (we have assumed that the 
upper band is completely empty) and 7k = ^{cos{kx) + 
cos{ky). In the limit W/t ^ 1 we have itk — >■ 1, Wk — >■ 
-^1^ so that we have approximately 




Corrections to this will be of order J/W. 
Whereas the small-Fermi-surface theory is valid for elec- 
tron densities close to n = 1 electron/unit cell, the large- 
Fermi-surface theory is valid for n = 2 electron/unit cell. 
In the next step we will compute the ground state energy 
as a function of the electron density. 
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FIG. 13: The function a{5) defined in (|55p evaluated numeri- 
cally in a 4 X 4 cluster for different forms of the kinetic energy, 
i.e. for different values of the hopping integrals t, t' and t" . 
The values of the hopping integrals for the individual systems 
are given m Table [D The line is the function ((56)) . 



V. PHASE TRANSITION IN STRONG 
COUPLING 

In the preceding two sections we have given a strong- 
coupling description of the two phases with large and 
small Fermi surface. Our main goal thereby was to eluci- 
date how these phases can be characterized beyond sim- 
ple mean-field theory. Clearly it would now be desirable 
to compare the energies of the resulting ground states 
and discuss a possible phase transition between the two. 
It should be noted from the very beginning that this will 
necessarily involve some rather crude approximations be- 
cause the quasiparticle Hamiltonians for the two phases 
a formulated in terms of Fermions with a hard-core con- 
straint. The considerations in the following section thus 
will necessarily have a more qualitative character. 
We first consider the dominant term in the Hamiltonian, 
Hw and collect only terms (x W. In the phase with 
the small Fermi surface the number of quasiparticles is 
NS and since each quasiparticle contributes — (see 
Eq. EZl) we have {Hw) = -^EM. + o(VFO). The 'vac- 
uum' for the construction of the large Fermi surface state. 



Eq. (|39|) . contributes an energy of — 



3W 



From the 



quasiparticle dispersion (|49p we obtain a contribution of 
-N{1 + S)^+0{W°) because the total number of quasi- 
particles is A^(l -I- (5). And finally there is an additional 
constant term of 2N which comes from inverting the 
two a-Fermion operators in the first term on the r.h.s. 
of Eq. gZl) so that we obtain (Hw) = -^MK, To 
leading order in W the two phases thus are degenerate 
which is very different from the mean- field treatment. 
We thus consider the kinetic energy. In discussing the 
kinetic energy we encounter the problem that the phase 
transition necessarily occurs for a doping range where at 
least one of the two states with different Fermi surface 
volume is far from 'its' vacuum state (these correspond 
to (5 = for the small Fermi surface state and (5 = 1 for 
the large Fermi surface state), so that relaxing the re- 
spective hard-core constraint between the quasiparticles 
will no longer be justified. To obtain at least qualitative 



results, we proceed as follows: we expect that for a finite 
density of quasiparticles the hard-core constraint leads to 
a reduction of the total kinetic energy. For a single-band 
Hubbard model we can study this reduction by exact di- 
agonalization of a small cluster. More precisely, defining 
the ground state energy of a single-band Hubbard model 
with Coulomb repulsion U and n electrons by Eq^\u) 
we can compute the function 



ain/N) - ^0 



oo 



(55) 



where N for the time being denotes the number of sites 
in the cluster. This gives the reduction of the kinetic 
energy due to the hard-core constraint. By evaluating 
both energies in the same 4x4 cluster we may expect 
to cancel out the shell-effects which inevitably dominate 
the kinetic energy of a finite cluster. In the numerical 
computation of a{n/N) we have moreover imposed the 
additional restriction to use only states with total spin 
S = 0. 

We have performed this procedure for different forms of 
the kinetic energy i.e. for different values of the hopping 
integrals t' and t" to (1, l)-like and (2, 0)-like neighbors. 
Figure [T2] then shows the resulting a{S) for different com- 



Ht 


t 


t' 


t" 


1 


1.00 


0.50 


0.40 


2 


1.00 


-0.50 


0.40 


3 


1.00 


0.50 


-0.40 


4 


1.00 


-0.50 


-0.40 


5 


1.00 


2.00 


2.00 


6 


1.00 


-2.00 


2.00 


7 


1.00 


2.00 


-2.00 


8 


1.00 


-2.00 


-2.00 



TABLE I: The values of the hopping integrals for the calcu- 
lation of the renormalization of the kinetic energy in Figure 
1131 The first column gives the number by which the dataset 
is labelled in Figure [131 
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binations of t' and t" . The respective values of t' and t" 
for each combination are given in Table HI As can be 
seen from Figure[l3]Q!((5) is relatively independent on the 
precise form of the kinetic energy, i.e. it seems to be an 
almost miivcrsal function of the particle density. a{5) 
can be fitted quite well by a simple quadratic function 
which depends on a single parameter 

a{5) = 1 + A(5 - (1 + \)5^ (56) 

where 5 — n/N . The fit gives the value A = —0.2. 
We now assume that the same reduction factor remains 
valid also for the more complicated Hamiltonians pro- 
duced by the above strong-coupling theories. It should 
be noted that this cannot be completely wrong because 
the values a(0) = 1 and q;(1) = are known and the func- 
tion may be expected to be slowly varying near 5 = and 
rapidly varying near 5=1. 

We therefore evaluate the kinetic energy for different elec- 
tron densities of the quasiparticles in the states with large 
and small Fermi surface by computing the kinetic energy 
in the absence of the hard-core constraint and then cor- 
recting by the factor a{5). Thereby for the large- Fermi- 
surface phase 5 is not the physical density of the conduc- 
tion holes, but the density of quasiparticles 

k,(T 

Let us stress that it is clear from the very beginning 
that this procedure must lead to a level crossing for a 
small value 6c of the physical hole concentration: in the 
absence of the correction factor a{5) the large-Fermi- 
surface-phase always has a lower kinetic energy because 
the renormalization of the hopping integrals is much 
weaker in this phase, compare Eqs. (|47l) and (I^Hl) . On the 
other hand, for (5 — >■ 0, the density of the quasiparticles 
(5 « 1 — (5 for the large- Fermi-surface phase approaches 1 , 
so that the correction factor a{5) approaches zero, wheras 
the density of quasiparticles for the small-Fermi-surface- 
phase is 5 = (5 and a(5) is close to 1. This gives a level 
crossing even if only the kinetic energy is taken into ac- 
count and we moreover expect this level crossing to occur 
for small 5. Figure [14] which shows 5c obtained by nu- 
merical calculation as a function of W/t confirms this. 5c 
decreases with increasing W/t which can be understood 
as a consequence of the decrease of the bandwidth in the 
small Fermi surface state with increasing W/t^ compare 
(EZl). 

Lastly, the d — d exchange term Hj should be considered 
but since we obtain a phase transition already from the 
kinetic energy and since the exchange energy is only a 
small correction for physical parameters we restrict our- 
selves to a qualitative discussion. For the small-Fermi 
surface phase we have the contribution from the exchange 
energy of the localized spins amongst each other, which 
was 2(1 — 2i5)Jxio- The discussion in section III sug- 
gested that xio is close to the value for the Heisenberg 
antiferromagnet, which is a lower bound for the nearest 
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FIG. 14: Hole concentration 5c where the phase transition 
from small to large Fermi surface occurs as a function of W/t. 
The small Fermi surface is realized for 5 < 5c. 



neighbor spin correlation function of localized electrons. 
The contribution from exchange between the quasiparti- 
cles, see Eq. ([M)) is cx 5"^ and since the kinetic energy 
favours a level crossing at small 5 - at least for larger 
W/t - this contribution is not important. For the large 
Fermi surface phase we first have a contribution from the 
virtual pair creation of triplets, Eq. ([5T|) . Since this is of 
order J"^ /W and will be suppressed because the density 
of quasiparticles is close to 1 this is not important either. 
Second, we have the contribution ([55)1 or ((M)) from the 
exchange between the quasiparticles themselves which is 
cx J. Since there the impact of the hard core constraint 
is hard to estimate we cannot definitely say which phase 
is favoured by the Heisenberg exchange - however, this is 
a small correction anyway. 

All in all we thus expect a phase transition between 
the states with small and large Fermi surface also in 
the strong coupling description. With the approxima- 
tions outlined above this transition should occur for any 
value of J/ i, which is different from the mean-field theory 
where a minimum value of J/t was required. Of course, 
quantitative agreement between mean-field and strong 
coupling theory may hardly be expected. 



VI. CONCLUSION 

In summary, the doping induced transition between 
phases with large and small Fermi surface in the 2D 
Kondo lattice model augmented by a Heisenberg ex- 
change between the localized spins has been studied by 
mean-field theory and in a strong coupling theory. Mean- 
field theory produces a phase diagram which has a rough 
similarity with that of cuprate superconductors: for low 
doping and low temperature the localized spins do not 
contribute to the Fermi surface volume but form a de- 
coupled spin liquid with pronounced singlet pairing. The 
spin liquid corresponds to a bond-related order parame- 
ter A(id with d3^2_y2 symmetry which describes singlet 
pairing between localized spins on nearest neighbors. 
This phase likely corresponds to the pseudogap phase 
in the cuprates. In mean-field theory the Fermi surface 
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of the conduction holes is a pocket around T which is 
unphysical due to the absence of any couphng to the lo- 
cahzed spins. One might conjecture, however, that the 
couphng to the antiferromagnetic fluctuations of the spin 
liquid would create hole pockets near (±-|,±-|) (see e.g. 
Refs. [6l|,|63) as possibly observed in ARPES[1[3. At 
higher doping and low temperature the localized spins do 
contribute to the Fermi surface volume and thus create a 
heavy Fermion-like state with a large Fermi surface and 
an enhanced band mass. This corresponds to the ovder- 
doped regime in the cuprates and is associated with a 
complex on-site order parameter Acd, which describes 
coherent local pairing between conduction electrons and 
localized spins. 

At intermediate doping and low temperature there is a 
phase where both order parameters coexist. This phase 
appears to be a dx2-y2 superconductor and has a large 
Fermi surface with a dj;2_y2 gap. Superconductivity 
emerges because the singlet pairing between the localized 
spins is transferred to the mobile conduction holes by the 
coherent Kondo pairing. All in all the phase diagram thus 
shows a certain analogy with that of the cuprates. 
In mean-field theory the superconducting transition is 
very different in the underdoped and overdoped regime: 
whereas in the overdoped regime we have a conventional 
2nd order transition with a dj.2_j^2-like gap opening on 
the large Fermi surface, the transition on the under- 
doped side corresponds to the emergence of a gapped 
large Fermi surface whereas the Fermi surface takes the 
form of hole pockets above the transition. This may actu- 
ally be consistent with experiment [1^. Mean-field theory 
moreover finds the transition to be 1"* order on the un- 
derdoped side of the superconducting dome. This is on 
one hand not really expected for a superconducting tran- 
sition but on the other hand an example of a I''* order 
transition involving competing order parameters [l^. If 
the transition really were first order an interesting possi- 
bility emerges: namely if the surface energy between the 
two degenerate phases were negative - as is the case in 
type-II-superconductors - this would imply a tendency to 
form inhomogeneous states in the underdoped region j53j. 
In fact one peculiar feature of underdoped cuprates is 
their 'granularity' [STt - lsgj . 

One experimental feature which is not reproduced by the 
present theory is the opening of a gap in an apparent 
large Fermi surface at the pseudogap temperature as ob- 
served by Hashimoto et al.^^. In the high-temperature 
phase the present mean-field theory predicts a complete 
decoupling of localized and conduction holes. 
To further elucidate the nature of the states with differ- 
ent Fermi surface volume and the existence of a phase 
transition we have also performed a strong-coupling cal- 
culation. In this theory, operators which create the eigen- 
states of a single cell from a suitably chosen 'vacuum 
state' are treated as effective Fermions. The advantage of 
this procedure is that the Kondo coupling, which should 
be the largest energy scale in the problem, is treated es- 
sentially exactly due to the 'pre-diagonalization' of a sin- 



gle cell. This theory then also yields a phase transition 
between states with different Fermi surface volume. As 
opposed to the mean-field theory, the localized and con- 
duction holes are approximately coupled to a singlet in 
both phases, so that the expectation value of the Kondo 
exchange is the same for both phases. In the phase with 
small Fermi surface realized at low doping the Kondo 
singlets - or Zhang-Rice singlets in the case of cuprates 
- form the quasiparticles so that the phase of a given 
Zhang- Rice singlet is determined by its momentum k and 
it is these momenta which form the small Fermi surface. 
In the large Fermi surface phase the Kondo singlets have 
a uniform phase and the momenta which form the Fermi 
surface are carried by sites with spin 1/2 and either 3 
or 1 holes (whereby the density of sites with 3 holes is 
oc {W/t)^). 

In both, the mean-field and the strong coupling theory, 
there occurs a transition between the two states with 
different Fermi surface for low carrier concentration. An- 
other complication which adds to the experimental com- 
plexity and which has been addressed not at all in the 
present paper, is a nematic ordering in the spin liquid 
plus hole pocket phase. 
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VII. APPENDIX A 

In this Appendix we want to address various issues 
related to the representation of the local singlets and 
triplets in terms of spin-1/2 Fermion operators. We note 
first that for example the states aj^^j^'o) and aj^^l'I'o) are 
orthogonal because their scalar product is proportional 
to (*o|^+|*o) =0. 

Deviations from ideal Fermion behaviour appears in the 
overlap of states with two Fermions. For example we 
have 

+ + 4 
t t 4 

4(*o|a»,t^jM,T^\i,t4tl^o) = ^ + 

4(*o|a»,t^M,i^\i,;aItl*o) = 1 - 

-+ + 4 
4(*o|a»,T^-,2,T^,2,t^ltl*o) = ^"3^»J 
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The factors of 4 on the left hand side thereby correspond 
to the prefactors in equation (IM)) . These relations would 
be consistent with those for ideal Fermi operators if 
the spin correlation function Xij = 0. If Xij is short 
ranged, the local singlets and triplets thus behave like 
Fermion operators for 'most distances'. We expect 
that the same will hold true for states with more than 
two Fermions provided that they are pairwise more 



distant than the spin correlation length C. This will 
be a reasonable assumption in the limit of low density 
that we are interested in. Furthermore, the neglected 
overlaps - being 'four particle overlaps' - would create 
an interaction between the quasiparticles rather than 
changing their dispersion. 
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